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We present a detailed description of biopolymer translocation through a nanopore in the presence 
of a solvent, using an innovative multi-scale methodology which treats the biopolymer at the micro- 
scopic scale as combined with a self-consistent mesoscopic description for the solvent fluid dynamics. 
We report evidence for quantized current blockade depending on the folding configuration and of- 
fer detailed information on the role of hydrodynamic correlations in speeding-up the translocation 
process. 



Biopolymer translocation through nanoscale pores 
holds the promise of efficient and improved sensing for 
many applications in biotechnology, and possibly ultra- 
fast DNA sequencing [3, 0, Q . Recent advances in fab- 
rication of solid-state nanopores 0, || have spurred de- 
tailed experimental studies of the translocation process, 
with DNA as the prototypical biopolymer of interest Q . 
Computer simulations that can account for the complex- 
ity of the biomolccule motion as it undergoes translo- 
cation, as well as its interaction with the environment 
(the nanopore and the solvent), are crucial in elucidat- 
ing current experiments @, Q an( i possibly inspiring new 
ones. Here, we study the dynamical, statistical and syn- 
ergistic features of the translocation process of a biopoly- 
mer through a nanopore by a multiscale method based 
on molecular dynamics for the biopolymer motion and 
mesoscopic lattice Boltzmann dynamics for the solvent. 
We report evidence for quantized current blockade de- 
pending on the folding configuration (single- or multi-file 
translocation) in good agreement with recent experimen- 
tal observations [7[. Our simulations show the signifi- 
cance of hydrodynamic correlations in speeding-up the 
translocation process. 

Nanopores are an essential element of cells and mem- 
branes, controlling the passage of molecules and regulat- 
ing many biological processes such as viral infection by 
phages and inter-bacterial DNA transduction Q. The 
last two decades have witnessed the emergence of artifi- 
cial solid-state nanopores as potential devices for sensing 
biomolecules through novel means [f|. One of the most 
intriguing possibilities is ultra-fast sequencing of DNA 
by measuring the electronic signal as the biomolccule 
translocates through a nanopore decorated with elec- 
trodes Q. While this goal still remains elusive, a num- 
ber of detailed studies on DNA translocation through 
nanopores have been reported recently [Tj, [8( . These ex- 
periments typically measure the blockade of the ion cur- 
rent through the nanopore during the time it takes the 



molecule to translocate, which provides statistical infor- 
mation about the biomolccule motion during the process. 

Numerical simulation of the translocation process pro- 
vides a wealth of information complementary to experi- 
ments, but is hindered by the very large number of par- 
ticles involved in the full process: these include all the 
atoms that constitute the biomolecule, the molecules and 
ions that constitute the solvent, and the atoms that are 
part of the solid membrane in the nanopore region. The 
spatial and temporal extent of the full system on atomic 
scales is far beyond what can be handled by direct com- 
putational methods without introducing major approxi- 
mations. Some universal features of translocation have 
been analyzed by means of suitably simplified statisti- 
cal schemes [Toj . and non- hydrodynamic coarse-grained 
or microscopic models [H QJ, El or other mesoscopic 
approaches [14j. Many atomic degrees of freedom, and 
especially those of the solvent and the membrane wall, 
are uninteresting from the biological point of view. The 
problem naturally calls for a multi-scale computational 
approach that can elucidate the interesting experimental 
measurements while coarse-graining the less important 
degrees of freedom. 

We have developed a multiscale method for treating 
the dynamics of biopolymer translocation [l5j and per- 
formed an extensive set of numerical simulations, com- 
bining constrained molecular dynamics (MD) for the 
polymer motion with a Lattice-Boltzmann (LB) treat- 
ment of the solvent hydrodynamics . The biopoly- 
mer transits through a nanopore under the effect of a 
localized electric field applied across the pore, mimicking 
the experimental setup [8(. The simulations provide di- 
rect computational evidence of quantized current block- 
ade and confirm the experimentally surmised multiple- 
file translocation: the molecule passes through the pore 
in a multi-stranded fold configuration when the pore is 
sufficiently wide. The simulations offer detailed informa- 
tion about several experimentally difficult issues, in par- 
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ticular the role of hydrodynamic correlations in speeding- 
up the translocation process. 

A three-dimensional box of size N x h x N y h x N z h 
lattice units, with h = Ax the spacing between lattice 
points, contains the solvent and the polymer. We take 
N x = 2N y , N y — N z ; a separating wall is located in 
the mid-section of the x direction, x = hN x /2. We use 
N x = 100 and N = 400, where N is the total number 
of beads in the polymer. At the center of the separating 
wall, a cylindrical hole of length Ihoie = 10/i and diameter 
dp is opened. Three different pore sizes (d„ = 5h, 9h, 17 h) 
have been used in the current simulations. Transloca- 
tion is induced by a constant electric field acting along 
the x direction and confined to a cylindrical channel of 
the same size as the hole, and length l p = 12h along the 
streamwise (x) direction. All parameters are measured in 
units of the lattice Boltzmann time step and spacing, At 
and Ax, respectively, which are both set equal to 1. The 
MD time step is five times smaller than At. The pulling 
force associated with the electric field in the experiments 
is q e E = 0.02 and the temperature is ksT/m = 10~ 4 . 
The monomers interact through a Lennard- Jones 6-12 
potential with parameters a = 1.8, and e = xlO -4 and 
the bond length among the beads is set at b = 1.2. The 
solvent is set at a density Plb = 1, with a kinematic 
viscosity vlb = 0.1 and a drag coefficient 7 = 0.1. 

We chose the separation d between the beads to be 
equal to the persistence length of double-stranded DNA, 
that is 50 nm, and define the lattice spacing to be 
d/1.2 = 40 nm. The hole diameters is 3 Ax. The re- 
pulsive interaction between the beads and the wall (with 
parameter a w = 1.5 Ax [l?]]) leaves an effective hole of 
size equal to ~ 5 nm. Having set the value of Ax, we 
choose the time step so that the kinematic viscosity is 
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expressed as: v w = vlb^t, with v w the viscosity of 
water (10 -6 m 2 /s) and vlb the numerical value of the 
viscosity in LB units; this procedure gives At ~ 160 ps, 
with vlb = 0.1. In order to ensure numerical stability, 
the relation ^At < 1 must be satisfied. Having estab- 
lished the value of At, we need to adjust the value of the 
drag coefficient accordingly, 7 < 6 • 10 9 sec _1 . This is sig- 
nificantly smaller than an estimate of the friction based 
on Stoke's law for DNA [l8| , which is equivalent to an un- 
dcrdamped system, or an artificially inflated bead mass. 
This approach is consistent with the coarse graining of 
the time evolution in the coupled LB-MD scheme. 

We focus on the fast translocation regime, in which the 
translocation time, t x , is much smaller than the Zimm 
time, which is the typical relaxation time of the poly- 
mer towards its native (minimum energy, maximum en- 
tropy) configuration. This corresponds to the strong- 
field condition q e Eb/kT > 1. In this regime, simple one- 
dimensional Brownian models [l^ | , or Fokker-Planck rep- 
resentations, cannot apply because the various monomer 
units do not have time to de-correlate before completing 
translocation. The ensemble of simulations is generated 
by different realizations of the initial polymer configura- 
tion, to account for the statistical nature of the process. 
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FIG. 1: Number of resident beads with time for three dif- 
ferent pore sizes d p = 5h, 9h, 17h (h = lattice spacing) and 
TVo = 400. (a) The fastest (minimum time, blue), slow- 
est (maximum time, red) and average speed (most probable 
time, green) translocation events; the insets show the cur- 
rent blockade for the duration of an event with average speed 
(green) with the current normalized to the open pore value 
(1). (b) Histogram P(N r ,t) of the distribution of JV r with 
time: short-time trajectories show multi-file character, reach- 
ing up to q ~ 2 — 8 in the initial stage of the translocation 
depending on the pore size; long-time trajectories show little 
departure from the single-file configuration. 



Initially, the polymer is generated by a three-dimensional 
random walk algorithm with different random numbers 
for each polymer configuration and one bead chosen ran- 
domly constrained at the pore entrance. Then, the poly- 
mer is allowed to relax for ~ 10 4 molecular dynamics 
steps without including the fluid solvent in the relax- 
ation, while keeping the bead at the pore entrance fixed. 
We define as time zero (t = 0), the time after the relax- 
ation, when the fluid motion is also added, the pulling 
force begins to act and the translocation process is ini- 
tiated; at this moment the bead at the pore entrance is 
also allowed to move. At this stage we do not include any 
electrostatic interactions within our model for reasons of 
computational simplicity. As far as the biopolymer mo- 
tion in the bulk of the solvent is concerned, this may 
actually be a good approximation of experimental condi- 
tions with high salt concentration, which leads to strong 
screening of electrostatic interactions. The situation at 
the pore region may require more refined treatment, be- 
yond the scope of the present work. 

In Fig. [TJa) we present the number of pore-resident 
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beads N r (i) as a function of time, for a narrow (d p = 5h), 
mid-sized (d p = 9h) and large (d p = 17 h) pore, with h 
the mesh spacing of the lattice Boltzmann simulation, for 
representative (fastest, slowest and average speed) tra- 
jectories. Simulations are repeated over an ensemble of 
400 realizations of different initial conditions and for to- 
tal polymer lengths up to No = 400. Time is measured 
in units of tf , the time it would take for the polymer to 
translocate if the monomers were to proceed in single- 
file configuration at the drift speed; this speed is given 
by ve = q e E/^m, with q e and m the charge and mass 
of the monomer, E the external electric field and 7 the 
hydrodynamic drag. This gives if = bN /vE = 12iV 
and the number of monomers in the pore for single file 
translocation is N± = 10 for the parameters used here. 

Fig. [3a) clearly shows the highly non-linear dynam- 
ics of the translocation process: In the initial stage of 
the translocation, the nanopore gets populated, with the 
number of resident monomers significantly overshooting 
the single-file value N\, the horizontal dashed lines at 
heights qNi indicating q-file (q = 1,2,3,...) transloca- 
tion. The range of q explored by the translocation tra- 
jectories grows approximately with the cross-section of 
the pore, going from q ~ 2 for the smallest pore d p = 5h 
up to q ~ 8 for the largest one d p = 17ft-. Note that 
these values correspond to about half the maximum al- 
lowed q-number, q max ~ d p /b. The fastest events corre- 
spond to the largest q value observed, while the slowest 
events correspond to essentially q = 1 throughout the 
translocation. It is also noteworthy that the transloca- 
tion time typically exceeds the single-file value, tf , except 
for the fastest events; for the most probable events q ~ 2 
for all pore sizes indicating that conservative monomer- 
monomer interactions produce an effective slow down 
compared to a single Langevin particle subject to a con- 
stant electric drive and frictional drag 7. 

Fig. nja) also presents the current blockade in all three 
pores for the most probable event in each case, which is 
the event with a translocation time close to the peak of 
the distribution over all translocation times. The current 
blockade is proportional to the number of monomers in 
the pore per unit area and appears to occur in well de- 
fined steps (quantized). Specifically, these blockades are 
calculated from the difference between the area of the 
resident beads, n(a/2) 2 , and the total area of the pore, 
7r(d p /2) 2 . In order to investigate the quantization of the 
current blockade we monitored the distribution of N r (t) 
at various time frame intervals of 100 steps. The res- 
ident monomers block the current across the channel, 
so that N r (t) conveys a direct measure of the current 
drop associated with the biopolymer passage through the 
nanopore. The corresponding histograms P(N r ,t) for 
three pore sizes are shown in Fig. [IJb). At early times, 
these histograms exhibit a multi-peaked structure, which 
is a clear signature of multi-file translocation. As time 
passes, the multiple peaks recede in favor of a single-peak 
distribution, close to the single-file value Ni = 10. This 
was found to be a stable-attractor for every simulated 



50 




1 1 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 

translocation time t x (t 1 units) 

FIG. 2: Scatterplot of the average resident number N r ver- 
sus translocation time (in units of tf) for the ensemble of 
translocation events for three values of the pore diameter, 
dp = 5h, 9h, 17h and N = 400. 



trajectory, indicating that the tail of the polymer always 
translocates as a single-file. 

Collecting all results for the average number of resident 
monomers N r as a function of the translocation time t x , 
for the three pore sizes studied, we find a simple rela- 
tionship, shown in Fig. [5J Experiments Q have reported 
that the average number of resident atoms N r in each 
translocation event varies approximately inversely with 
the duration of translocation t x : 




N r (t) dt = N r t x oc N = const. (1) 



The single-file asymptote Ni = 10 (q = 1) at long- 
times, t 3> tf , is evident. The short-time asymptote, 
reaching up to 4 < q < 5, corresponds to ultrafast 
translocations (t < tf) occurring in the case of the 
large-diameter pore, d p — 17h. These results arc in- 
tuitively reasonable, since large resident numbers im- 
ply that more monomers cross the pore per unit time, 
hence the translocation becomes faster. The results also 
support the notion of N r (t) as a measure of the time- 
rate of the translocation, dNx/dt oc N r , from which the 
inverse-proportionality between N r and t x is a direct con- 
sequence of J Q tx [dNT/dt] dt = No — const. In this expres- 
sion, Nx{t) is the number of translocated monomers at 
time t. 

The simulations reveal that solvent correlated motion 
makes a substantial contribution to the translocation en- 
ergetics. The role of hydrodynamic correlations is best 
highlighted by computing the work done by the moving 
fluid on the polymer (we call this the synergy, Wh) over 
the entire translocation process as compared to the case 
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FIG. 3: Statistical distribution of the work performed by the 
hydrodynamic and electric field during translocation events. 
The vertical dotted line corresponds to the work Wf done by 
the electric field on polymers that translocate single-file. 

of a passive fluid at rest: 

rU N ° 

W H {t x )=l y^Viit) ■ Vi{t) dt (2) 
J o i=l 

where Vi is the velocity of monomer i and Ui is the fluid 
velocity at the position of monomer i. For the sake of 
comparison, it is also instructive to contrast Wh with 
the corresponding work done by the electric field 

W E (t x )=q e Y^E-Vi(t)dt (3) 
J ° <=i 

where the sum extends over the resident monomers only, 
since the electric field is applied at the pore region only. 
These statistically averaged values of Wh and We reveal 
a number of interesting features (see Fig. [3]). First, Wh is 
always positive, clearly showing that hydrodynamic cor- 
relations provide a cooperative background, as compared 
to the case of a passive "ether" medium (u = 0). Second, 
we observe that the W E has a much narrower distribution 
of values than Wh, reflecting the ordered structure of the 
biopolymer as it passes through the nanopore, as com- 
pared to its off-pore morphology. It is useful to introduce 




FIG. 4: Left panel: a typical two-folded polymer configura- 
tion (d p — 9h), at time where 65% of the No = 400 total beads 
have already translocated from right to left; colored contours 
show the magnitude of the corresponding hydrodynamic syn- 
ergy field (only five of the nine wall-layers are shown). Right 
panel: a single-file translocation event for a narrow and shal- 
low pore (dp = 3h, l p — lh) with 60% of the beads translo- 
cated and the corresponding magnitude of the synergy. 



the work done by the electric field on molecules which 
translocate single-file and proceed through the pore at 
speed v E , Wf = q e ENiv E tf = bq e ENiN . In the ab- 
sence of any other interaction, a g-file translocation at 
speed ve would complete in a time t x (q) = tf/q un- 
der an electric work qWf . In the present simulations, 
Wf = 0.12A , thereby Wf = 48 for N = 400. Inter- 
estingly, the distribution of We values is highly peaked 
at a value very close to Wf. The observation that 
W E ~ Wf implies that qv x (q)t x (q) ~ vstf and since 
the simulations show that t x (q) > tf , the conclusion is 
that v x (q) < VE/q, indicating that collective motion of 
the monomers slows down the process. 

A major asset of numerical simulations for the study 
of translocation processes is the direct access to visu- 
alization of the morphology of the translocating chain. 
As an example, we show in Fig. 5] a typical "snapshot" 
at a time when about 65% of the monomers have al- 
ready passed through the pore of a translocating 2-foldcd 
chain of Ao = 400 beads. In the same figure we show for 
comparison an event for the same length, but for single- 
file translocation (unfolded chain) through a very narrow 
(dp = 3h) and shallow (l p = 1) pore. In addition to the 
polymer conformation, we show isocontours of the mag- 
nitude of the hydrodynamic synergy density 

w ff (r;<)=7 M*) • ( 4 ) 

ieB(r) 

which is a local (in both space and time) version of the 
total synergy Wh defined in Eq.Q, with B(r) a grid 
cell centered around location r = (x,y,z). The con- 
tours of Wn{f) illustrate the cooperative nature of the 
hydrodynamic field, with regions of high co-moving flow 
surrounding the translocating polymer and assisting its 
motion. This is suggestive of the notion of an "effec- 
tive" polymer, dressed with the hydrodynamic synergy 
field, which acts as a self-consistent lubricant, helping 
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FIG. 5: Average translocation time as a function of polymer 
length, with (closed circles) and without (open squares) hy- 
drodynamic interactions. Colors correspond to different pore 
diameters d p — 5h (green solid line), 9h (red dotted line), and 
17h (blue dashed line). Black triangles indicate the value of 
the single-file translocation time tf for each value of No with 
hydrodynamics (all numbers are scaled to the value of tf for 
N = 400). 



the polymer to negotiate a faster passage through the 
nanoporc. 

We further investigate this issue by inspecting the av- 
erage (over an ensemble of 400 realizations) translocation 
time, (t x ), as a function of the polymer length, with and 
without hydrodynamics. The results are shown in Fig. [5] 
hydrodynamics consistently accelerates the translocation 
by roughly 30% percent. More intuitively, hydrodynam- 
ics literally re-normalizes the diameter of the pore: as 
is clearly visible in Fig. O a pore of diameter d p = hh 
for a bare polymer (without the hydrodynamic field) is 
essentially equivalent to a pore of almost double diame- 
ter dp = 9h for the hydrodynamically-dressed polymer. 
In order to assess the degree of correlation between the 
translocation dynamics of the "dressed" polymer versus 
the actual one, we have measured the translocated spe- 
cific synergy (synergy per monomer), defined as 
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N T (t) 



N T (t) 



Ui(t) 



(5) 



with Nx(t) the number of translocated monomers at time 
t. Clearly, any implicit time-dependent functional de- 
pendence of the form {wnit)) = (wn(NT(t))) would in- 
dicate that mass and synergy translocate in a synchro- 
nized manner. We find that the ratio (u ■ v)/kT ~ 5, 
reflecting the fact that the solvent locally "follows" the 
monomer and providing a measure of the relative impor- 
tance of synergistic versus thermal forces. Our results 
show that (wait)) is essentially constant throughout the 
translocation process. This implies a direct proportion- 
ality between the translocated synergy and the number 
of translocated beads and supports the notion that the 
"dressed" and the actual polymer proceed in full syn- 
chronization across the nanopore. 

In conclusion, by using a new multiscale methodol- 
ogy based on the direct coupling of constrained molec- 
ular dynamics for the solute biopolymers with a lattice 
Boltzmann treatment of solvent dynamics, we have been 
able to confirm a number of experimental observations, 
such as a direct relation between quantized current block- 
ades and multi-folded polymer conformations during the 
translocation process. In particular, the simulations re- 
veal an intimate connection between polymer and hy- 
drodynamic motion which promotes a cooperative back- 
ground for the translocating molecule, thus resulting in a 
significant acceleration of the translocation process. Such 
an acceleration can also be interpreted as the outcome of 
a renormalization of the actual polymer geometry into 
an effective one, more conducive to translocation. This 
opens up exciting prospects for the development of op- 
timized nano-hydrodynamic devices based on the fine- 
tuning of hydrodynamic correlations. As an example, one 
may envisage multi-translocation chips, whereby multi- 
ple molecules would translocate in parallel across mem- 
branes with an array of pores. The optimization of such 
devices will require control of solvent-mediated molecule- 
molecule interactions to minimize destructive interfer- 
ence between translocation events. 
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